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ABSTRACT 
Forecasts  of  ARIMA  processes  are  generally  made  using  the  Difference 
Equation  form.   This  is  the  approach  favoured  by  Pox  and  Jenkins  and  most 
subsequent  authors.   The  purpose  of  this  note  is  to  emphasise  that  ! 
Integrated  For::;  of  the  forecast  enjoys  some  important  advantages  derived  from 
its  explicit  use  of  the  Eventual  Forecast  Function  (EFF).   A  brief  review 
the  procedures  for  obtaining  all  the  necessary  components  of  the  Integrated 
Form  is  given,  and  a  new  and  lirect  method  is  derived  for  evaluating  the  extr? 


t  ime  s . 


In  their  book,  Box  and  Jenkins  (1970)  strongly  recommend  that  forecasts 
of  ARIivlA  processes  be  made  using  the  Difference  Equation  fern  because  it  is 
the  simplest  approach.   All  subsequent  textbooks  have  endorsed  this  view  to 
the  extent  that  very  few  of  them  consider  forecasting  in  the  Integrated  Form 
in  any  detail.   A  notable  exception  to  this  is  the  recent  book  by  Abraham  and 
Ledolter  (1983)  which  contains  sere  useful  detailed  discussion  of  the  role  of 
the  Eventual  Forecast  Function  (E7F)  in  generating  forecasts.  The  purpose  :f 
this  note  is  to  indicate  that  the  Integrated  Form  does  have  value  and  is  worthy 
of  consideration.   V/e  begin  with  a  brief  review  of  the  derivation  of  the 
necessary  components  for  this  :  :r~   of  *xo  forecast. 

T'HE  I^JTE jRAmZD  FOEd.1  OF  '~q:r  FOPfZCAST 

Suppose  {."■'  }  is  an  Ail  I  MA  rrocess  satisfying  a(3)X  =  8(5^ 

t  t        h 


where  B  is  the  Eackshift  inerator  which  is  defined  b'. 


A .   -  A  . 


p    ...       q 

where  a(3)  =     a,  3""  .  0(E)  =   V   8,  3~  ,  a,  =  p~  =  1   and  all  seasonal  and 
k=0  k=0 

difference  terms  ere  included  in  the  i—     an  i   8-  operators.   For  lead  tio.es 

F 

T  >  q  ,   y   a,,X+(T-k)  =  0  ,  where     X  (i)   is  the  forecast  made  at  time  t  :f 

k=0  " 

X     and  takes  the  value   X   .   if   i  <  0  .   The  solution  cf  this  differenc 
t+i  t+i        - 

equation  is  the  EFE.   It  car.  be  written  as 

? 
X  (T)  =  I     b^f,  (T)  (i; 

.-r =  - 

where      f    (T)      are    deterministic   functions    of      T      and   may    include    polynomials, 
k 

exponentials,    sinusoids   and  products   of    these.      They   ,oay   also   be    dummy 
variables   generating    a    seasonal   pattern. 

As   a   representation   of    the   T-steo   ahead   forecast    (1)    above    is    '/slid   '    v 
cniv      T   >   :,!   =   o    -    z    .      If      M  >    0    .    then   for    lead    times      T    -    I,.. ,\1 


b.':\  (T)    -  cL    .e,    .  (2) 


is  the  sequence  of  one-step  ahead  forecast  errors,  i.e. 

et  -  xt  -  h-lil]  ■ 

Box  and  Jenkins  (1970 )    show  that  the  current  value  of  b,  =  (b, ,b0,...b  ) 


may  be  obtained  from  b  1      via  a  linear  equation  b_  =  Lb    +  h e ,  .   The 

t  —  -i-  X      X —  -L       X 


matrix  L  effects  the  changes  in  the  coefficients  in  revising  the  time  origin 

from  (t-1)  to  t  ,  and  can  be  obtained  as  L  =  R.  F,1+1  ,  where  F>,  is  the  (pxp) 

_]_ 
matrix  with  (i,j)th  element   f.(M+i)  .   The  vector  h  =  F,,  (k.  ,  where 

i  k 

■  '  =    (oM  -,  ,-,._.  ,...,■:,   )   and  cl>v  is  the  coefficient  cf  B 

r\3)    =  j3(B)/a(B)  ,  the  usual  moving -average  r   -   ntation  of  the  ARIMA 

process. 

V/e  r.\  -   n   t      t    at   the   revision  matrix     1     is   block   diagonal.      Each  block 

corre;        ids   to  a   real    (repeated)    linear  factor   or   conjugate  pairs   of   complex 

factors    in  a(B)    .      foe   effect    of   this    structure    is    that   individual  components 

liter      .    re vis  -  pen  of   the  others.      Thus,    for 

example,    any   linear         ad  a     rial   fad  :rs    can  be   obtained  at  each  time     t 

•i  :-:■:.   only   their   values   at   time    (t-1         .  I    the   latest  forecast   error     e_,_    . 

Deta:  the   derivation   of    the   revi    .'     .   equati  'or   the   linear  trend  and 

;onally    differenc  Lei  forecasts  a  von   in  McKenzie    (19S4b). 

]  ration   and   inversion  of   the  matrix     F   . 

If    ,ve  wish   tc  'ely    in.'     '  components   of   the   SFF  but  the 

or  .,    we   must   use    (1)    or    (2)    above.      If 

brief!;       '        :  and  .  enkins    \  n  .      )pears   to  be 
a   messy  proth  .  fever,    v/e  pr       >r.t  here   a     ;ene  -al  solution   relat- 

?f  i  ■ '     .   .■     .  •  -'    the   model. 


The  derivation  of  this  result  is  given  in  an  Appendix  at  the  end  of  the 

There  it  is  shown  that  cL   =  d        ,  (j  =  1,2,...,M-T  ;  T  =  1,2,...,U). 

Thus,  it  is  necessary  to  obtain  only    {cL  n  :  T  =  1,2,...,M}.  and  these 

1  ,U 

easily   derived  from  the   following  matri 


a        a     .      a     _    . 
p       p-1       p-2 


a  a     ,    . 

P  P-l 


1 


£ 


i,o 


3,o 


L1,0 


" 


p+1 


p+3 


o: 


The  triangular  nature  of  these    .  tions  evidenl  L;    kes  then;  particular!; 

to  solve.   This  is  especially  true  when  ■    note  that  ".  =   q  -  p   is  very  rar  Ly 

large. 

Ecx   and  Jenkins   give    some   non-s  xamples       .  ir  hook   and   a 

sonal  one   is   given   in  McKenzie    (19  expressing    th  -  ;  .   .        tirline   model 

forecasts    in      irrent    levs    ,     gradient    and   seasonal   factor  form.      Vie   give   one 

other  brief   example  here   to   illustrate    the   evalua    Ion  of  the     1-c   efficie 
Consider 


1    -    :b)(l    -   B)X     =    (1 
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;      -(!+*] 


0  : 


"i.  : 

*-*■      s 

-. 

■ 

Ld2,0j 

: 

i.e. 


-  ,    \        -; 


'1, 


The  Integrated  Form  of  the  forecast  is  now 


Xt(l)  =  bx  +  b2?  +  dljQet  +  d^0et_1  , 


.  ,  t      u  d 

X^{2  I     =  0-,  +  0.,  i       +  c 


^;0    o 


XJT)  =  bf  +  b 

0  1 


6     V  >  9 


where 


b^ 

1  0 
0   6 


fb*-1 


,**l 


J  U. 


(l-8)/(l-9)   le. 


(e-<r)/r(l-<i> 


nrc  ""'"'COT  "V ! 


,  .  .  rS       U  i  J 


*  p  r»  o  in 


i:ue  —u  l^cl  i,j.u;j 


,;ie 


sinniest  mathenatica 


•-SJ  >   - 


is  important  to  note  tnat  it  need,  net,  be  the  mos' 


eiiicient  computationally. 


r  3   ^-»o 


:fficiency  of  the  two  forms  depends 


or:  wna" 


-A_q   sro 


i ore casting   nas   already    ~een 


i  ;t    and   how   often.      This    aspect   of   ARLMA 
;ed  by    jodciphin    fiQr"t)    and  more    recently  by 


icr.ersie    ;  iVc^a  i  . 


he    twe    forms    correspond    to    two    distinct    approaches 
Difference    Ec nation   procedures    generate    for  ^casts    re cursive It   f 


\^ ^   *  -    -        — .   .- -  •—  q .-   «._,..  ..  ^    .  . .         .  -     . .    ~        r  orm     or    t        ot r er   riar  d 

involves    the    calculation   cf    a   set    :f    values   which   can  be   used    to   construct 

As   an   illustration,     consider   the  ARL'.-AK  ,.  ,2      proc        .      zqy.  and   Jenkins 
( 1970,pi4-6- )    discuss    forecasting    this   process    in    detail   using   both   approaches, 
so  we    simply    rei   '         :e    the    firms    here.      The    Integrated   Form   is    considered   first, 
It   has    r.vc    distinct   carts: 


-  ,    -_-.v 


'.;■',-   -jT\: 


,-I 


=     0  + 


'  >  \ 


(ii)  Forecast  Equation 

Xt(k)  =  o1q   +  kbjj  k  >   1  (5) 

t     t 
There   is   a  fixed  computational   investment   in  generating  b_  =    (b    /c,  ) 

before  any   forecast   can  be   obtained.      Once   this    is   achieved,  ver , 

forecasts   for  any   lead  times   are   easily    ierived  using  equation    (5).      All  the 

storage   requirements   relate  to    the    re  n  equations    (4),    and  he        '     ralve 

For   the  Difference   Equation   form,    the   forecasts    ars  ited   directly: 

Xt(l)    =     Xw(2)    + 

x,  (2)  =    ;<         )    -  .  ^  (o) 

X.  (k)    =  2X, (k-1)    -  1    ;k-2)  k  >   3  (7) 

to  u 

Note  that  there  are  no  i    visi    squat ioj  ;  ,v  ic   a  not  generate  for      . 

However,  it  is  clear  that  the  for'     s  of  lead  times  k=l,2  and  .:  play  a  role 

in  equations  (6)  and  (7J  similar  to  that  :>f  b,  in   luations  (4.1  and  (5?. 

Another  important  point  here  is  that  to  obtain,  forecasts  for  any  lea:  "' 

using  (6)  and  (7)  we  must  first  generate  the  forecasts  for  all  shorter  lead 

times. 

To  il   trate  the  differences  1       a   t'A   app]  .aches  .ve  consider  some 

simple  fc     ;t  ■  ;  scenarios.   Supj       aish  to  forecast  for  lead  tia.es  k  = 

1,2,...,T.   The  cc   .   tional  requir    ts  are  clearly  comparable. 

integrates  Form  requires  less  storage  to  generate  t     recasts,  but  the 

ion   •  laticns  necessitate  a  little  r:.ore  arithmetic  caan  the  use  ::'  (&) 

and  (7).   Suppose  new  that  the  lead  tia.es  of  interest  are  not  ;on    Ltive. 

For  example,  u  .   ,vish  forecasts  for  lead  times  k=  ,2,3,  .  .'-      ■' 

Using  thf  integrated  Form,  ;  -     It        te  only  t;      Lx  forec 


has  been  obtained.   However,  the  Difference  Equation  form  requires  the 
generation  of  all  18  forecasts.   Alternately,  we  may  routinely  forecast  for 
lead  times  1,2  and  3  and  occasionally  require  forecasts  for  others,  e.g.  6,12, 
18  etc.   Again,  it  is  easier  tc  generate  these  via  (4), (5)  than  (6), (7). 

Another  prediction  of  common  interest  is  the  cumulative  forecast,  i.e. 

the  forecast  of  Y,  (T)  =     X.   .   It  is  given  by  Y  (T)  -  '   X  (k).   Using  the 

t        ,_1u  +  K  o        ,   ,  t 

K  -  j.  £■  - 1 

Difference  Equation  form  we  must  generate  X.  (k)  for  k=l,2,...,T  .and  then  sum 
them.   For  the  Integrated  Form,  we  can  use  a  single  forecast:- 

m 

y\  (T)  =   7  ibt  -  kb?)  -   Tb+0  +  ^T(T+l)b? 
>  =  I  °  °    - 

These  comments  apply  equally  -.veil  to  other  situations  in  which  linear 

functions  of  future  values  are  to  be  predicted.   They  are  usually  more  easily 

treated  _  ■  h  -  the  Integrated  Form. 

There  is  no  s  :    tion  here  that  one  act  reach  is  always  better  than  the 

-  - '-■  g  r>   ~c  —  '"■■=•■"  ~ "" ~  ~    pa  nr    has  certai1^  ^dvant^^es  in  certain  situations.   The 

investment  it.  the  Revision  Fcuations  c:    the  Integrated  Form  yields  ^r0^31 


;«  ucj 


generation   and  a   variety   of    forecasting   problems   are 
thus   more   efficiently    treated  by    this    ai     .         a.      To   assess   the   cost   cf    such   an 
investment,    art-?    that    the    linens:    a    >f    the   vect  >r   b     is   the   order   oi     z'r.e 

•ressive-dif f erence  r;  tor    in   the   ARIMA  process,    e.g.    (p+d)    in  an 

ARIi/A(p,d,q]    process.      Thus,    revision    :f   b  +    is    roughly    2omparable    to   the 
production   of    forecasts    for   the    first    (p+d)    lead   times    using   the   Difference 
Equation    form.       Clearly,    if   no    leaa   tia.es    beyond    (p+d)    are    to   be   predicted, 

iffer   ■  Equation   form   is   more   efficient    in  general.      However,    when  lead 

t  ime  s    an    exes     '      f    :  a ■*■  : .    d re    rec uire d  i tter   smc-roach    ie rends      as   we   ha  -~  ■ 


On  the  other  hand,  there  is  a  most  important  advantage  enjoyed  by   the 
Integrated  form.   It  lies  in  the  area  of  interpretation.   This  is  an  aspect  of 
forecasting  which  should  not  be  underestimated.   V/e  .-ray  recall  Stem's 
comments  (1974)  that  managers  who  require  forecasts  may  be  prepared  to  accept 
trends  and  seasonal  effects  because  these  correspond  to  familiar  ideas.   This 
view  is  merely  one  of  many  such  cris  de  coeur  still  heard  by  forecasters.   The 
true  value  of   a  forecast  is  not,  invested  solely  in  its  accuracy  but  also  in 
its  credibility.   The  former  can  be  assessed  only  after  the  event  forecast, 
whereas  the  latter  will  let ermine  whether  the  forecast  is  used  at  all.   The 
modelling  of  time-series  as  ARIMA  processes  is  now  a  common  practice,  thanks 
to  the  availability  of  a  variety  :f  powerful  computer  packages.   !;everthei- 
a  major  hurdle  for  most  lsers  is  still  the  practical  interpretation  of  the 
forecasts. 

The  reason  that  the  Integrated  form  enjoys  an  advantage  here  is  that  i" 
generates  the  forecast  via  the  Eventual  Forecast  Function  'IF?].      This  is  a 
linear  combination  of  deterministic  functions  of  the  lead  time  T.   For  fixe  : 

thp  fnnor  i  ,-ir  o  qT*°  Vn  ~\\m    j/offlv   hnt  thp1'  r  cop^f  ^  c\  ert  ^  "n  'hoe   '"  'a1-1 
combination  a  cart  with  each  new  observation.   Thus,  In  the  AEIMA(  ),2,2) 
examclp   the  IFF  is  -^iven  bv  ecuaticr  (5)  .   The  ie termini stic  functions  are 

f1(T)=i  and  fp(T)=T.  Their  coef f icien+s ,  the  components  of  bf,  are  revised 
with  each  new  observation  via  equations  (4).  The  new  observations  allow  the 
forecast  to  adapt,  while  preserving  the  basic  structure  :f  the  EFF. 

This  structure  is  most  important  because  +ue  individual  ietermir :'  stic 
functions  whin,  appear  in  the  IFF  can  o^t^n   te  readily  interpreted  in  terms 
of  familiar  concepts  such  is  trends,  growth  and  seasonality.   It  is 
interesting  to  note  that  it  is  the  practice  of  the  •  .        id               :recastij 
systems  such  as  Holt -Winters  or  Exponential  f    thing  *   model  the  :'  t  as 


jus:,  such  a  linear  combination  of  deterministic  functions.   Surely  ,  the 
success  and  -.vide  acceptability  of  these  systems  is  at  least  partly  due  to  the 
fact  that  their  predictions  are  cast  in  a  familiar  form.   In  addition,  more 
recently  popularized  systems  such  as  the  Eayesian  models  of  Harrison  and 
Stephens  (1976),  or  the  Kalman  Filtering  approach  discussed  by  Harvey  (1983), 
are  also  constructed  around  these  same  familiar  components  of  trend  and 
seasonality,  etc.   On  the  coder  hand,  none  of  the  software  available  for 
modelling  and  forecasting  ARL'.IA  processes  can  generate  forecasts  in  this 
readily  understood  form.   Nevertheless,  such  interpretations  are  available  vi- 


■  -     '  v-    _j   »^    - 


'. 


jt    tne   or '    m   toe   AK1..1A  '.  ere  cast 
f-xt^aci  '       nonit  3rec      and   cro^ected. 

As   an   illustration,    .ve   consider    :he  monthly   sales    of   US   houses    (January, 
lco:;-~-ceooer;  1975 )  .      Tne    iata   are  given    in  Abraham  and  Ledolter    (1933),    and 
tpp   fitted   mo  * ~'j ^    *-"^r'^'    i 

:   this   model    .vere    obtained   from  the    iota   using   backforecasting.      The 
components    cf    the   Integrated  Forra  of    the   forecast  were   evaluated  and   the 


(1-2    d-5"  )X^  =    ( 1-0.23 )( 1-0.833"    )a. 


err 


r-  p  r*  o  q    :  p  f;      >-»o  T 


ry    revision  eauaticns   obtamec    isin^    ::.      ;  r    secures     lescrioea   in 


::.\  ;:^.  a- 


•-q  • 
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!,"      +    0.74 
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t  _- 
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■      .0935e 


:"|       ?  1  1 

i.    ,  _    ,    .    .    .     ,   -Li 


(9) 


Y     ■  >  .      =    '-   J 


(10 


.7.  .>;..     ,\  , 


rv  1  e   =  X 


X.    ^(1)    evaluated,    the 


revision  equations  yield  current  values  for  the  non-seasonal  component  bx  = 

— t 

(b    ,bn  ),    and  the   seasonal  comuonent  S,    =    (3. ,S0 , . . . ,3, _ ) .      These   cor 
o      1  — t  i '    2  12 

have   useful  and  readily  understood   interpretations.      In   the  non-seasonal 

t  t 

component,    b     is    the   current   level  of    the  rrecess,    and  b.    the   current    gradient 

'   o  i  "J 

of  the  linear  trend,  i.e.  the  predicted  rate  of  change  in  level  per  month.   In 

the  seasonal  component,  S  is  the  additive  seasonal  factor  :~-r   month  (t+>).   It 
predicts  the  amount  by  which  the  data  will  deviate  from  the  process  level  '.-: 
months  from  now  'time  t). 

Exhibit  ]  .'  out  here  . 
Exhibit  1  iispiays  the  lata  for  the  years  1968-75,  and  the  corresponding 
values  of  the  current  level.   The  values  of  the  gradient  art  not  displayed, 
follow  a  path  similar  in  shape  to  that  o:    the  level,  though  with  a  iifferent 
vertical  scale.   In  January  1968,  with  a  lata  value  of  15,  we  find  the  current 
level  is  4.1.8  and  the  gradient   .C67.   Thus,  in  this  month,  the  non-seasonal 
component  c:    the  forecast  for  T  months  ahead  is  given  by  (4-1.3  -  C.067T). 
Clearly,  as  the  data  evolved,  tnis  proved  to  be  optimistic,  as  the  crane  during 
the  next  two  years  appears  to  be  iewn wards.   This  is  reflected  in  cue  January 
1970  gradient   .027,  vhich  is  smaller,  but  still  positive.   From  equation 
we  can  see  that  the  gradient  acaccs  only  slowly  with  each  new  ?vservaticn.  In 
January  1971  and  1973,  the  gradient  adapts  to  the  increasing  level,  with  valu 
0.173  and  ),236  respectively.   In  the  latter  month,  the  predicted  linear  trend 
in  the  process  T  mont]  ;  .head  is  (59.7  -  0.2367).   As  Exhibit  1  clearly  sh 
this  would  have  been  a  very  coor  lon£  term  sredicti   .   However,  by  January 
L974,  the  predicted  trend  is  (41. t  -  C.061T).   This,  in  the  event,  ray  be  a 
little     Lmistic.   It  is  Lmnortant  to  n    that  \  t  n.1        ave  we  been  a: 


quote  the  non-seasonal  component  without  reference  to  the  seasonal  one  but, 
from  equations  (8)  and  (9),  we  can  calculate  them  separately.   Thus,  for 
example,  we  could  forecast  using  the  Difference  Equation  form  and  monitor  the 
level  and  trend  directly  via  13). 

Exhibit  2  about  here 
The  values  of  the  seasonal  component  are  revised  with  every  observation 
and  can  be  plotted  in  the  same  way  as  the  non-seasonal  one.   Since  it  involves 
twelve  paths  and  they  are,  in  fact,  very  stable,  they  are  not  reproduced  here. 

It  is  just  as  revealing,  in  this  case,  to  plot  the  current  seasonal  component 

+ 
of  the  :ata,  i.e.  X  -  b  .   This  is  displayed  in  Exhibit  2.   There  is  a  very 

clear  ar  i  ~,crsister.~  seasonal  oattem.   The  months  March  until    rust  'ire 

always  high  and  approximately  the  same.   The  ether  months  are  equally 

consistent,  .vith  the  possible  exception  of  Eecember,  which  appears  to  be  on  a 

iownwara  trend.   It  is  interesting  to  note  that  when  a  subset  ;f  these  data  tup 

to  May  19751  was  analysed  -arlier  by  Killmer  and  Tiao  (1979]  those  authors 

chose  •.  "■-  .  -- 1  in  which  the  ....-,.   deterministic.   The 

applicability  ::'  such  a  model  is  reflected  in  Exhibit  -.   he  would  have  mud. 

more  confidence  In  prediction    f  the  .'-      L  com.!    r.z ,      _,  of  this  forecast 

than  the  n  n-         ne  , 

3U?.!MA?.Y 
The   purpose    ;f    this    paper    Ls    to   make   a    :ase    for    the   usefulness   of    the 
Integra  Form  i    'T-X.'.'f.   for     ;ast.  rief    review   of    its     ierivation   is 

presented.      .'.  new   method   for    Dbtaini;   ;  essential   part   of    it    is   given   in 

-  e  z  a  1 1 .  o  rm  '         *  ~~ '  •  !  c  ■ .    ~      '  *    ■  n   n  e  '■"  1 ':-  \  t  e  i    ir    t  hf3    lite  ra  t  u  r  e 

3Ut    it   ha;  'uJ     iharacteristies .  i  irive    from   its    direct    use    ::    the 

Eventual   Forecast  Function.      Ir.        rl  Lcuiar,    we   have  ised   two  points. 


11 


First,   it  need  not  be,  as  is  often  supposed,  computationally  inferior  to 
the  Difference  Equation  approach.   Secondly ,  it  offers,  via  the  2FF, 
considerable  advantages  for  interpretation  of  the  forecast.   In  bo        ,  the 
argument  is  illustrated  with  an  example. 
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APPENDIX 

Derivation  of  the  d-coef f icients 

First,  note  that  substituting  equation  (2)  into  the  usual  forecast 

revision  identity  X.  (T)  =  X,  ,  (T+I)  +  ilime,.  and  using  the  revision  equation 
J         z  t-i  it7         ° 

b^  =  Lb.  ,+  he  +  yields   d^   =  d      .  ,  (j  =  1,2,...,M-T  ;  T  =  1,2, ...,M). 

Thus,    it    is  necessary   to  obtain   only      {cl,        :    T  =   1,2,...,M}      and  this  ray  be 
achieved  as   folic '.vs. 

Applying   one  usual  conditional  expectation   arguments    to     X,      to  derive   T-step 
ahead  forecasts    for     T   =  p+l,    p+2,...,q   ,   yields 

?  .  '■-" 

7      a.  X    Ltc-r-k)     =  -  <r    =    1    0  M'  (11) 


-•—  '        L     ^ 


a,f.(k)    =  C,    the   left-hand  side   of    (11 


may   be  wri  ,ten  as 


x  ■  A,  ,  o+r-: 

Lj  f         t     * 


where  lj  =  rrax  .,  p+r-M  .   A  further  application  cf  equation  [Z'i   yields 


• 


x         :    ,  . e ,  .  (11) 

k   .-„  *r-k,i  t-i 


ivhi  :h  is  *       iated  to  the  right-  :  i  ;ide  ;f  (11 ) . 

Equating  :oeff icients  ::'  e   .   for  i  =  3,l,...,M-r  ,  ,ve  find 

t  -i 


(13) 


P+r— :Yi+i  for       i   <        <  M— r  ^   r    <  '.'.    .      From   these 

at  ions      1;      it    is   easy    t       ieri  atrix   eauati    .. 
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Exhibit    1.   US   House  sales,  Jan   58-Dec  75. 
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